Conductances, Conductance Fluctuations, and Level Statistics on the Surface of 

Multilayer Quantum Hall States 
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The transport properties on the two-dimensional surface of coupled multilayer heterostructures 
are studied in the integer quantum Hall states. We emphasize the criticality of the surface state 
and the phase coherent transport properties in the thermodynamic limit. A new, stable numerical 
algorithm for large scale conductance calculations in the transfer matrix approach is discussed in 
detail. It is then applied to a directed network model describing the quantum mechanical tunneling 
and impurity scattering of the multilayer edge states. We calculate the two-probe conductance in 
the direction parallel to the external magnetic held, its fluctuations and statistical distributions 
as a function of the interlayer tunneling strength. Using finite size scaling, the asymptotic scaling 
functions of the ensemble averaged conductance and the conductance fluctuations are calculated for 
a fixed aspect ratio and found to be in remarkable agreement with the analytical results obtained 
using the supersymmetric nonlinear cr-model approach. The conductance distribution is determined 
in the quasi-one-dimensional metallic, insulating, as well as the crossover regime where comparisons 
are made to that at the single-layer quantum Hall transition. We present, for the first time, a 
detailed study of the level statistics in the eigenvalue spectrum of the transfer matrix. Coexistence 
of metallic and insulating statistics is observed in the crossover regime, which is attributed to the 
emergence of a finite range level repulsion in the crossover regime, separating the metallic (Wigner- 
surmise) behavior at small level spacings from the insulating (uncorrelated or Poisson) behavior at 
large level spacings. 



PACS numbers: 73.20.Dx, 73.40. Hm, 73.33.-b, 72.15.Rn, 05.30.-d. 



I. INTRODUCTION 



When a two-dimensional electron gas (2DEG) is placed 
under a strong perpendicular magnetic field, it exhibits a 
remarkable set of low-temperature magnetotransport be- 
haviors known as the quantum Hall effect (QHE) 
The main part of the phenomenology can be summa- 
rized by the existence of (1) new stable phases of matter, 
i.e. the quantum Hall states, with vanishing dissipation 
and quantized Hall conductances; (2) continuous, zero 
temperature phase transitions between the quantum Hall 
states; and (3) extended, current-carrying, edge states in 
the quantum Hall liquids. These edge states are intrinsi- 
cally chiral due to the complete breaking of time reversal 
symmetry by the external magnetic field. 

Two-dimensionality has been regarded as a necessary 
condition for observing the QHE. Recently, the ques- 
tion concerning what happens to the physics, e.g. the 
three aspects mentioned above, associated with the QHE 
in materials having engineering dimensions between two 
and three has received substantial theoretical and exper- 
imental interests. Integer quantized Hall plateaus have 
been observed experimentally in weakly coupled 30 Q] 
and 200 || multilayer GaAs/AlGaAs graded structures. 
In the presence of finite interlayer tunnclings, the latter 
are the natural generalizations of the QHE above 2D. 
There are several important issues H^]: (1) the quan- 
tization condition for the Hall resistance; (2) the stable 
phases of matter and the nature of the phase transitions; 
(3) the new physics associated with the edge states phe- 
nomena in the quantum Hall states. 





FIG. 1. Schematics of a multilayer quantum Hall sam- 
ple. The dotted lines indicate the coupling between the 
edge states on the surface due to quantum tunneling. 

In this paper, we shall focus on (3), under the con- 
dition that the multilayer graded heterostructure is in 
the bulk quantum Hall state. In this case, the bulk lo- 
calization length is finite and very short. The edge 
state in each layer thus decouples from the bulk. New 
physics arises because the finite interlayer tunnelings on 
the sample boundaries couple the edge states supported 
by individual layers together to form an interesting 2D 
system on the surface of the multilayer Hf/J. This is 
shown schematically in Fig. |j} At this stage, one assumes 
that the electron-electron interactions can be ignored and 
the edge states can be treated as independent electrons 
undergoing impurity scattering and interlayer tunneling 
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as described in Fig. |J. 

The three most noteworthy novelties of the surface 
state under such settings are as follows, (i) It is a 2D chi- 
ral system because of the unidirectional edge states, (ii) 
As such, the backscattering mechanism of 2D weak local- 
ization is effectively suppressed. The surface electrons are 
in a critical state, in the sense that the localization length 
£ is linearly proportional to the linear dimension of the 
surface [^). (iii) The transport is intrinsically anisotropic. 
Let us consider such a sample surface as shown in Fig. [I], 
with dimensions L in the interlayer direction and a layer 
circumference C . In the presence of disorder, the chiral 
nature of the edge states keeps the in-plane transport bal- 
listic with a velocity vb- The typical time required for 
an electron to circumnavigate the sample is thus given 
by a ballistic time tb — C/vb- The transport along the 
direction parallel to the magnetic field (z-direction) be- 
comes, on the other hand, diffusive characterized by a 
diffusion constant D which is controlled by the strength 
of the interlayer tunneling t. This leads to a diffusion 
time tjj = L 2 / D across the sample in the z-direction. 



1D Insulator 



1 D Metal 
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FIG. 2. Schematics of the crossover between different 
phase coherent transport regimes on the surface (see text). 
The dashed line corresponds to a constant circumference 
C, its intersections with the solid crossover lines determine 
Lq and £ respectively. 

The transport property in the magnetic field direction 
(z) has been a recent subject of considerable theoreti- 
cal |6|j7|,|^-[l^] and experimental interests jL5|. The most 
striking feature is the emergence of several regimes of 
phase coherent transport in samples where the phase co- 
herent length 3> L, C. These regions are most clearly 
exhibited in the plane spanned by L and the surface as- 
pect ratio, 



(1) 



as shown schematically in Fig. [|. (1) The localization 
length of the critical state, £ ~ C, separates the metallic 
regime for L < £ from the quasi- ID insulating regime 



for L £. This crossover is shown by the horizontal 
line at A — const, in Fig. 0, where L = £. (2) Next, 
equating the ballistic time tb and the diffusion time td 
leads to a new length scale Lq ~ \[C. Associated with 
the latter is another crossover between a chiral 2D and 
the quasi- ID metallic regime, as shown in Fig. [2] by the 
curved line with A oc 1/L on which L = Lq. In the 
2D metal regime, the electrons typically diffuse across L 
before ballisticly traverse across the circumference C . (3) 
Finally, the regime of small L and large A is expected to 
be dominated by ballistic transport. 

For fixed dimensions L and C, and thus fixed aspect 
ratio A, the conductance in the z-direction is expected to 
show smooth crossovers between the three regimes as a 
function of the interlayer tunneling t. Alternatively, for 
a fixed t, the crossovers take place with varying sample 
geometry. The dashed line in Fig. ^| shows an example 
of the latter, along which C is kept fixed. As L increases 
and sweeps through Lq and £, transport properties cross 
over from those of a 2D chiral metal to those of the quasi- 
1D metal and finally to those of a quasi-lD insulator. 

In this work, we are interested in the transport proper- 
ties along the field direction in the thermodynamic limit, 
i.e. when C — > 00 and L — » 00. A large scale numeri- 
cal transfer matrix algorithm will be discussed and ap- 
plied to calculate the two-terminal conductance and its 
mesoscopic fluctuations in phase coherent samples. The 
behaviors in the thermodynamic limit will be obtained 
via finite size scaling. Note that the transport behaviors 
of the surface state in the thermodynamic limit depend 
crucially on how the latter is taken. In particular, since 
£ oc C, the thermodynamic state reached along a path 
with a finite but fixed aspect ratio A = const, corre- 
sponds, for any interlayer tunnelings, to the ID metallic 
or insulating regimes, where L/Lq ~ L/yC — > 00, as 
shown in Fig. ||. The 2D chiral metallic regime, on the 
other hand, can only be reached in the thermodynamic 
limit with a vanishing aspect ratio (see Fig. [2]) A ~ L~ a 
and a > 1. The scaling forms for the ensemble averaged 
conductance and its variance have been written down by 
Gruzberg, Read, and Sachdev (GRS) @H, 



e*(L L 

<ff> " r U% 



and similarly, the 2n-th central moment follows, 



(2) 



(3) 



Here, the crossover lengths L and £ depend on C and 
the strength of the interlayer tunneling t. Note that the 
scaling variable L/£, with £ being the localization length 
in a finite system, contains the aspect ratio (A) depen- 
dence in the scaling limit where £ oc C p6[ . 

In the quasi-lD regime, it was pointed out by GRS 
that the scaling functions T and T2 should follow those 
in a quasi-lD wire jl3 14 1. The asymptotic forms of 



2 



the latter have been obtained exactly by Mirlin et al. 
using the supersymmetric nonlinear cr-model fl7j| . The 
exact scaling functions in the chiral 2D regime are not 
known. Recently, they have been studied perturbatively, 
and the leading order corrections to the scaling function 
have been obtained |]l4| . The first comparison of numer- 
ical calculations with these analytical results was made 
by Cho, Balents, and Fisher j28|. In addition to study 
the behaviors of the mean conductance and its variance, 
we will also analyze the behaviors of the conductance 
distribution function and the level statistics in different 
regimes and near the crossover. 

The rest of the paper is organized as follows. In section 

II, we will describe a directed network model for trans- 
port on the surface. We will discuss the properties of the 
transfer matrix for chiral transport and compare it to the 
one used in the 2D quantum Hall transport. In section 

III, we will discuss a new numerical algorithm for large 
scale conductance calculations. The method for study- 
ing the level statistics is also explained. We will then 
present our numerical results in section IV and conclude 
in section V together with a discussion of the open issues. 

II. DIRECTED NETWORK MODEL 

To facilitate a calculation of the transport properties 
in the presence of a smooth-varying impurity potential, 
we model the motion of the electrons on the surface (see 
Fig. £l) by the directed network (DN) model shown in 
Fig. |j introduced by Chalker and Dohmen [|6|JT^] . The 
propagation of the chiral edge states is represented by 
the directed links in the DN. The impurity scattering 
is accounted for by letting the wavefunctions accumulate 
random Aharonov-Bohm phases along the links |lj]] . The 
interlayer tunneling between the edge states take place at 
the nodes of the DN. 




i 
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FIG. 3. The directed network model. The links with 
directed arrows indicate the propagation of the edge state 
in individual layers, whereas the interlayer tunneling takes 
place at the nodes. The magnetic field is along the 
z-direction. 

We define the dimensions of the DN in Fig. ^ as fol- 
lows: L along the horizontal (z) direction, corresponding 
to L columns of nodes where tunneling between adja- 
cent layers takes place among the L coupled layers; C 
in the transverse (x) direction, which amounts to a layer 
circumference of C/2 rows of nodes. For simplicity, the 
lattice constants of the DN will be set to unity such that 
L and C are dimensionless numbers. Periodic bound- 
ary conditions are applied in the transverse direction, 
whereas two semi-infinite ideal leads will be attached to 
the ends for the calculation of the two-terminal conduc- 
tance in the z-direction. 

The difference of the DN from the Chalker-Coddington 
network model for the integer quantum Hall transport 
|]l9| is that the probability current on all the links is uni- 
directional. As a result, the structure of the transfer 
matrix describing a quantum tunneling event at a node 
is different. As shown in Fig. ||, at each node, there 
are two incoming and two outgoing modes. Let us de- 
note the associated probability amplitudes by ^L,in = ^2, 
^L,out = Zi, Zii. la = Z 3 , and Z R . ou t = Z\. The transfer 
matrix must be defined in a way such that the transfer 
process is multiplicative throughout the DN, 

/ -^L.out \ _ rp I -^R.out \ (*\ 

\ Z h . ln ) ~ \ Z KM ) ' [V 

where T is a 2 x 2 transfer matrix. The current 
conservation at each node implies that T must leave 
J2i=L,R.\ z iAn\ 2 invariant, i.e., 

|^L,in| 2 — |^L,out| 2 = |^R,out| 2 — |^R,in| 2 - 

Thus, T has noncompact unitary symmetry, T S {7(1, 1). 
With a choice of gauge, one has, 

rp ( sinh# cosh#\ 

^ cosh 6 sinh(9 ) ' 1 ' ' 

controlled by a real tunneling parameter 9. Notice that 
detT = —1, i.e. for the chiral surface transport, T be- 
longs to the noncompact unitary group with determinant 
negative unity. In general, is determined by the Fermi 
energy and the local potential. We will take to be the 
same on all the nodes, corresponding to uniform inter- 
layer tunnelings. We found that a randomly distributed 
9 about a finite mean value across the DN is irrelevant for 
the asymptotic transport properties of the critical surface 
states. This situation is reminiscent of the one encoun- 
tered in the ordinary network model description of the 
integer quantum Hall transition fl9[[ . In the latter case, 
the randomness in the tunneling amplitude at the saddle 
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points of the random potential was found to be irrelevant |2(j,^lJ . 

The scattering S'-matrix, on the other hand, transforms the incoming amplitudes to the outgoing ones at a node, 



r t 

t -r 



J R,in 



where t and r are the transmission and the reflection coefficients respectively. From Eqs (||) and (||), one finds, 

1 



t = 



cosh 9 ' 



= tanh#. 



(6) 



(7) 



In the rest of the paper, we will study the transport properties as a function of the single-node interlayer tunneling 
strength i 2 . 

It is instructive to compare the transfer matrix in Eq. ([bJ) to the one describing the quantum Hall transport within 
the layer @, 



_ _ / cosh^ sinh# 
QH ~ I sinh6> coshfl 



(8) 



The transfer matrix in this case belongs to the noncompact unitary group J7(l, I) with determinant det Tqh = 1. Thus, 
the two problems are very different and belong to different ensembles of random matrices, except in the pathological 
limit of 9 — > oo, i.e. r — > I and t — > 0, where T — Tqh- 

Using Eq. (0) as the building block, one can construct the C x C transfer matrix for the i-th column of the DN in 



Fig. § 



T t = 



V 



o 









J4>c 



T 



(9) 



Here fa S [0, 2ir) are the random U(l) phases on the links of the DN. The transfer matrix for the i ± I-th column has 
a slightly different structure, 



/ e^i sinh6> 



Ti±i — 





e #3 



cosh 9 \ 



T 



Ke^'c coshfl 











T 



(10) 



e^'c sinh 9/ 



for periodic boundary conditions. Here fa € [0, 2ir) are again the random U(l) phases. 

across the entire DN implies that T must leave 



The total transfer matrix for the entire DN is given by 
the matrix product, 



(ii) 



a/2 

i=i 



C/2 



I^L,out| 2 ) — E (I^R.outl 2 _ I^R,in| 2 ) 



Correspondingly, the scattering S'-matrix has the form 



Thus T forms a noncompact unitary group U(C/2, C/2) 
that has the following metric preserving properties: 



r t 
t -r 



(12) 



T f ET = — E, 
TET f = — E, 



(13) 
(14) 



where t and r are the (C/2) x (C/2) transmission and 
reflection matrices respectively. 

Now we discuss the symmetry properties of the transfer 
matrix product in Eq. (11). The current conservation 



where the C x C matrix E can be written as a direct 
product ((g)), 



E = r 2 



1, 
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of the 2x2 Pauli matrix r z and the ((7/2) x (CI 2) identity 
matrix 1. Notice the minus sign in Eqs (|13| ) and (14). 
This is a consequence of defining the transfer matrix that 
is multiplicative across the DN, as in the single-node case 
in Eq. (^). As a result, detT = — 1, i.e. the noncompact 
unitary group is not volume-preserving. 

What we will show next is that the relevant hermitian 
matrix T^T entering the conductance formula (see sec- 
tion III) has the desired properties. Substituting Eq. (H~ 
to the left hand side of Eq. (|l3|) for E leads to, 



(15) 



(T t T)£(T t T) = S. 



Thus, T^T forms the volume-preserving U(C,C) group. 
Denoting the eigenvalues of the latter by {A^}, i = 
1, . . . , C, in descending order, Aj > A i+1 , we next show 
that the eigenvalues form reciprocal pairs, 

A, • Xc^+i = 1. (16) 

Under a unitary transformation, 

P = U^T^LI, U+U = 1, 

the matrix P has the same set of eigenvalues {A^}. Fol- 
lowing Ref. |1, Eq. © implies that 



ptjp = j } j = U f SU. 
With the following choice of U, 

U = -j= (r ® 1 - ir x <g> 1) , 
the metric preserved by P in Eq. (Jl7|) becomes 



(17) 



/ 1 

-1 



J = ir v 



\ 



V 



/ 



Thus the P matrices belong to the symplectic group 
Sp(C). As a consequence, the eigenvalues of P, and thus 
those of T^T, form reciprocal pairs as in Eq. (|T6|). This 
property will be used in the conductance calculations in 
the following section. 



III. NUMERICAL ALGORITHM 

To determine the two-terminal conductance in the 
magnetic field direction (z), two semi- infinite ideal leads 
are attached at the left and right ends of the disordered 
DN in Fig. ||. The conductance is then given by the Lan- 
daucr formula [E3lE4|, 



g=-Tr[t% 



(18) 



where t is the (C/2) x ((7/2) transmission matrix defined 
in Eq. (|l2|). Using the unitarity of S and the symplec- 
tic property of T^T, it is straightforward to show that 
Eq. ( |18|) can be written in terms of the transfer matrices 

il, 



g = X Tr 



TtT+(TtT) L +2-l 



(19) 



Further simplication can be made by writing the ordered 
eigenvalues ({Xi}) of T^T as, 



A, 



-2 7 , = 



,,C/2, 



i =£7/2+1,. 



(20) 



Evaluating the trace in Eq. (|l9|) leads to, 



C/2 



-Y 



\ cosh (7j) 



(21) 



In the remainder of the paper, we will set the conductance 
unit e 2 /h = 1. 

The two-terminal conductance can be calculated from 
Eq. (|2l|) or Eq. (|l^) for a given sample with a disorder re- 
alization. The conductance distribution is then obtained 
by repeating the calculation for a large ensemble of sam- 
ples with microscopically different disorder realizations. 
The ensemble averaged conductance, conductance fluc- 
tuations, and in fact all statistical properties of the con- 
ductance can be obtained from the conductance distri- 
bution. Nevertheless, this conceptually straightforward 
procedure has been hindered by the numerical difficulties 
involved in evaluating the product of a large number of 
random transfer matrices (see e.g. Eq. (|ll|)). Unlike the 
localization length, which can be obtained in the quasi- 
1D limit of a long stripe or a cylinder using Oseledec's 
theorem [^6| , the conductance of a C x L sample requires 
the knowledge of the product matrix of 2L random trans- 
fer matrices of dimension C x C (Eq. (|l9| )), or equiva- 
lently, that of all its eigenvalues (Eq. (|21[)) . For large L 
and C, direct calculation of the product is impossible as 
one quickly looses numerical stability in multiplying the 
transfer matrices due to the inevitable dynamical range 
problem. The latter refers to the inability of today's com- 
puters to keep within accuracy of exponentially growing 
and exponentially decreasing matrix elements. As a re- 
sult, the applicability of this method has been limited to 
small system sizes. 

Now we discuss a stable numerical algorithm that re- 
moves the above difficulty. We shall first prove a theorem, 
which is then applied to the evaluation of the eigenvalues 
7, in Eq. @. 

Consider a general hermitian matrix H with non- 
degenerate eigenvalues. 

Theorem 1 The n-th power of H, i.e. the matrix H n 
can be written in the form: 
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H" — U n D n R n , 



(22) 



where, in the limit of large n, (i) U n is a unitary matrix 
of which the columns converge to the eigenvectors of H ; 
(ii) D n is a diagonal matrix and the eigenvalues of H 
is given by D n D~^_ 1 in descending order; (Hi) R n is a 
right triangular matrix with unity on the diagonals that 
converges to a limiting matrix of that structure. 
Proof: Let us begin with the fact that the hermitian 
matrix H can be factorized according to, 



H = U1D1R1. 



(23) 



This UDR-decomposition is identical to the well-known 
Gram-Schmidt procedure for orthonormalizing the col- 
umn vectors in H. Apply H to Eq. (p3|), one has, 



H 2 = HUiDiRi = U2D' 2 R / 2 D 1 R 1 . 



(24) 



Here we have done another UDR-decomposition to the 
product HU\. Thus R' 2 is right triangular with unit di- 
agonal elements. Now we interchange R' 2 and D±, 



R' 2 D\ — D\R 2 , 



(25) 



where R 2 is right-triangular with non-vanishing ele- 
ments, 



Di{i,t) 

Thus Eq. (|24|) becomes, 

H 2 = U 2 D 2 R 2 , 

where 

D 2 = D' 2 D\, R 2 = R 2 Ri, 

are diagonal and right-triangular matrices respectively. 
Repeating the steps from Eq. ( p3| ) to Eq. (27) leads to 

H n = U n D n R n , 

and the following recursion relations, 

D n = D' n D n _ l7 

Rn — R'lRn-1 ■ 



(26) 



(27) 



(28) 



(29) 



(30) 
(31) 



Eq. (|29j) is a matrix-factorized form of the product ma- 
trix H n . The unitary matrix U n is necessarily well con- 
ditioned. We will show below that the right-triangular 
matrix R is well defined and converges to a limiting ma- 
trix of that structure exponentially fast with increasing 
n. The large variations in the size of the matrix ele- 
ments, which limit the dynamical range of the matrix 
multiplications, are contained solely in the elements of 
the diagonal matrix D n . In fact the Gram-Schmidt or- 
thonormalization procedure ensures that, in the limit of 
large n, the columns of the unitary matrix U n project 
onto the eigenvectors of H. 



Next, we demonstrate the convergence of R n and ob- 
tain the eigenvalues of H. The important point is that, in 
each step, the UDR-decomposition is done by orthonor- 
malizing the j-th column with respect to all columns of 
index i < j. As a result, the diagonal elements of D n 
are accumulated in descending order, i.e. D n (i,i) S> 
D n (j,j) for j > i. Carrying out one more multiplication, 



H 



n+l 



HU n D n R n — Un-i-iD'n^Rn^DnRn 



Interchanging R' n+ i with D n , 



R'n+lDn — D n R 



n+l) 



where 



K +1 (i,j) = R' n+1 (i,j) 



Dn(j,j) 

D n (i,i) ' 



(32) 



(33) 



(34) 



Since R(i,j) is nonzero for j > i, Eq. ( |34|) shows that the 
off diagonal elements (J > i) are suppressed exponen- 
tially with increasing n by the ratio of D n (j, j)/ D n (i, i), 
while the diagonal elements remain to be unity. Thus, 
when n is large, R'n+i approaches the identity matrix 
provided that there is no degenerate elements in D n . 
The convergence of the R matrix follows since R n +i = 
R" +1 R n — + R n for large n. Going back to Eq. (|3^), we 
have 



H n+1 — HU n D n R n — U n +iD n+ iR n 



Using the fact that R n +i — ► Rn, Eq. 

Ul +1 HU H = D n+1 D-' 



implies, 



(35) 



(36) 



Thus we have proven that for large n, the eigenvalues of 
H is given in descending order by the diagonal matrix 
£) n _|_i-D~ , and the columns of the unitary matrices U n 
and U n +i converge to the corresponding eigenvectors. In 
the absence of nearly degenerate eigenvalues, the proce- 
dure converges exponentially fast, such that in practice 
n need not be too large to obtain the desired accuracy. 

It is straightforward to construct an algorithm for eval- 
uating the eigenvalues of the transfer matrix product in 
the conductance Eq. ( ^l| ) following the steps of the proof 
described above. Specifically, we consider H as defined 
by 



H = T f T = TlTl ■ ■ ■ TjTi, ■ ■ ■ T 2 T X 



(37) 



The eigenvalues of H is obtained by the above UDR- 
decomposition of 



H n = 



(38) 



In practice, the UDR-decomposition need not be per- 
formed after every multiplication, as this is the most 
time-consuming part of the computation. It scales at 
best with C a , 2 < a < 3. The number of possible direct 
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multiplications without losing accuracy, Nd, depends on 
the size of the matrix elements in T^, thus the tunneling 
parameter t, and weakly on the size of the matrices. 

Algorithm: 

1. Carry out Nd steps of direct matrix multiplications 
in Eq. Q. 



2. Perform a UDR-decomposition as described above. 

3. Repeat the above procedure for a total of (n — 1) x 
(2L) multiplications for the convergence of the fac- 
torized matrices. 

4. For the last 2L multiplications, store and accumu- 
late the logarithm of the elements in the diago- 
nal matrix for each of the N u d r = 2L/Nd UDR- 
decompositions as 



2ji = logAj = ^2 lo S Di(i,i). 



(39) 



;=i 



For large enough n, {A^} gives the set of eigenvalues 
for H in descending order. 

5. To check convergence, apply H again to the unitary 
matrix at the end of the last UDR-decomposition. 
Repeating step 4, one has 



E 

i=i 



log D t 



(40) 



The difference between {7^} and {7,'} can be com- 
pared to the convergence criterion, and thus deter- 
mine the appropriate power n. 

6. The conductance can be readily calculated from 
Eq. @. 

We emphasize that the applicability of this stable algo- 
rithm for conductance calculations is very general. It is 
particularly useful for studying phase-coherent transport 
properties in disordered macroscopic and mesoscopic sys- 
tems in the thermodynamic limit. The only limitation in 
this case is the CPU time. It has been successfully ap- 
plied to study the conductance and conductance fluctua- 
tions in the 2D integer QHE [ plp^ . The specific choice 
of the two parameters, n and Nd, depends on the sys- 
tem under investigation (see section IV). In general, it is 
more convenient to choose Nd such that N u d r is an inte- 
ger. Otherwise, an additional UDR decomposition must 
be performed at the end of each set of the 2L multiplica- 
tions. As evidenced in the procedure leading to the proof 
of the above algorithm, essentially, only U and D need 
to be kept, and R can be discarded in each of the UDR 
decomposition. Thus, even the standard Gram-Schmidt 
procedure works for the present purposes. In the next 
section, we will apply this technique to study the sta- 
tistical properties of the transport on the surface of a 
layered quantum Hall state. 



IV. THE QUASI-1D REGIME 

To study the quantum interference effects in the con- 
ductance, it is convenient to introduce the Ohmic con- 
ductance in the classical limit, <? hm- The latter can be 
obtained analytically for the DN in Fig.|^ by summing 
over all possible non-crossing Feynman paths contribut- 
ing to the z-axis transport, while ignoring the nontrivial 
interference effects associated with the intersecting ones 
winding at least once around the circumference of the 
DN |§,Eq]. This leads to the averaged conductance, 



C 



9o = 



2Ll-t 2 (l-l/L) : 



(41) 



where t 2 is the interlayer tunneling amplitude defined 
in Eq. (Q). The L-dependence in the denominator in 
Eq. (^d]) is indicative of the ballistic contribution that 
dominates in samples with L <C t 2 /(l — t 2 ). Thus, the 
crossover from ballistic to diffusive transport is rather 
slow if t 2 is close to one. 

In the limit L — > 00, one recovers the Ohmic behavior 
in the z-direction from Eq. (f4l|), 



C 



1 t 2 



2 1 - t 2 ' 



(42) 



As a result, the dependence of measurable quantities 
on the microscopic parameter t 2 can be favorably re- 
placed by the dependence on the Ohmic conductance 
<?ohm or the conductivity a. For example, as discussed 
in Refs the 1D to 2D crossover length L in the 

scaling equations (^) and (|J) can be expressed as 



and the localization length £ as 



(43) 



(44) 



We consider square DNs with aspect ratio A = 1, un- 
less otherwise specified. The asymptotic behavior of the 
transport properties in the thermodynamic limit taken 
with L = C — > 00 will be determined using the finite size 
scaling (FSS) analysis. In this case, 



L 



= A\ 



(45) 



such that the system is always in the quasi-lD regime 
in the thermodynamic limit (see Fig. ||). In fact, this 
is true so long as the latter is taken with a fixed finite 
aspect ratio, A — const.. Nevertheless, 



L 
I 



A 
2^' 



(46) 



which allows one to study the crossover between the 
quasi- ID metallic (£/£ ^ 1) an d insulating (L/^ ^> 1) 
regimes as a function of the interlayer tunneling in the 
thermodynamic limit. 
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A. Conductances and Conductance Fluctuations 

For a given disorder realization, we calculate the con- 
ductance from the Landauer-formula in the form of 
Eq. (pl| ) using the algorithm described in the last sec- 
tion. Next the conductances of a large ensemble contain- 
ing 5, 000 to 10, 000 microscopically different samples are 
calculated. This is done for a sequence of sample sizes of 
8 x 8 to 96 x 96 to facilitate the FSS analysis. The typical 
parameters used are the following. The direct number of 
multiplications Nd ranges from 8 for large t 2 (approach- 
ing one from below) to 64 for small t 2 (approaching zero). 
The power n is typically less than 50 which is more than 
enough to suppress the systematic errors in the obtained 
eigenvalues below the statistical errors associated with 
ensemble averaging. The latter is typically less than 1% 
for our data. 

We first present the results for the ensemble averaged 
conductance. In the quasi-lD regime, the scaling func- 
tion in Eq. (^) becomes, in the asymptotic limit, 

< 9 >=r^,ooV (47) 

For fixed aspect ratio, Eq. ( f46| ) shows that L/£ is only a 
function of t 2 . 

As an example, the calculated < g(t, L) > is shown 
in Fig. | as a function of L at t 2 = 0.76 or £/£ = 0.31, 
which is in the metallic regime. The apparent size depen- 
dence is due to the corrections to scaling. To extract the 
asymptotic value, we perform a FSS analysis according 
to, 

< g(t, L) >=< g(t) > -Cirr-^J-, (48) 

where y; rr = 0.74 ± 0.03 is the dimension of the leading 
irrelevant operator that controls the corrections to scaling 
and £j rr is a non-universal constant introduced by the 
conjugate finite length scale. 




FIG. 4. Finite size scaling of the averaged conductance 
in the quasi-lD metallic regime at L/£ = 0.31. The solid 
line is a fit to Eq. @, with y irl = 0.74 ± 0.03 and the 
asymptotic value < g >= 1.60 ± 0.01. 

Following this procedure, the asymptotic values of the 
conductance are obtained for different values of t 2 . In 
Fig. pj, the conductance normalized by the Ohmic value, 
< g > /g hm, is plotted as a function of L/^ in accordance 
with Eq. (fl7]). The solid line is the analytical scaling 
function r(L/£) obtained by Mirlin et al. for a quasi- ID 
wire using the supersymmetric non-linear tr-model [ p~7[ | . 
The agreement throughout the entire regime is remark- 
able, which is in accordance with the analytical work of 
GSR |l3| . It is important to note that such an excellent 
agreement can only be obtained for the conductances in 
the thermodynamic limit by the FSS analysis of large 
sample sizes. In fact, the nature of the finite size correc- 
tions is quite different on the metallic and the insulating 
side of the crossover. Let us define 

f3(t) = -^<g(t,L)>. (49) 

One finds that (3(t) — > as L — > oo consistent with the 
fact that the system is critical for all t 2 ||. However, we 
find that there exists a well defined tunneling amplitude, 
t% ~ 0.37, such that (3(t) -> 0+ for t 2 > t%, whereas 
[3{t) — > 0~ for t 2 < t 2 r . It is thus natural to associate the 
existence of t 2 r with a sharp crossover from the metallic 
to the insulating regime. The corresponding crossover 
conductivity is a cl - ~ 0.3 and (L/£) cr ~ 1.7. These ob- 
servations provide an explanation of the systematic de- 
viations from the analytical curve observed numerically 
in relatively small systems p8[. 




2Ut, 

FIG. 5. The averaged conductance normalized by <7ohm 
as a function of 2L/£. The data points are the thermody- 
namic limit values obtained using finite size scaling. The 
solid line is the analytical curve for a quasi- ID wire ob- 
tained by Mirlin et al. 
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Next we present the results for conductance fluctua- 
tions about the ensemble-averaged values. Once again, 
we observe significant size dependence in the variance 
of the conductances. The latter is indicative of the 
fact that the universality of these conductance fluctua- 
tions is a property of the quantum transport at criti- 
cality pl|p7yi4] , and is outside the context of the con- 
ventional universal conductance fluctuations in diffusive 
metals p9[ | . Eq. (|^) has now the asymptotic form, 



< Sg 2 >= r 2 



, oo 



(50) 



The asymptotic values of the variance are plotted in 
Fig. ^| as a function of L/£ to display the scaling function 
r 2 in Eq. (p0|). The solid line is the analytical result of 
r 2 obtained by Mirlin et al. [[ItJ . As in the case of the 
ensemble-averaged conductance, the agreement for the 
asymptotic variance of the conductance is excellent for 
the entire quasi- ID regime. 
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FIG. 6. The variance of the conductance in the asymp- 
totic limit as a function of 2L/£. The solid line is the 
analytical curve of T2 for a quasi- ID wire obtained by 
Mirlin et al.. 

In the limit of small L/£, the asymptotic variance ap- 
proaches the universal value < Sg 2 >= 1/15, despite the 
vanishingly small fluctuations at small L where the trans- 
port is ballistic. In the opposite limit of large L/t;, the 
system is in the ID insulating regime where the absolute 
fluctuations decay exponentially. Nevertheless, it is im- 
portant to note that in the insulating regime, the averaged 
conductance itself decays exponentially |l7j , 



V< Sg 2 > 
<9> 



(52) 



deep in the quasi- ID insulating regime. 

We emphasize that this is the first time that con- 
ductances and conductance fluctuations are numerically 
evaluated in the thermodynamic limit in this quasi- ID 
regime. The agreement between the numerical data ob- 
tained from the DN model and the analytical results from 
the supersymmctric nonlinear a-model is, as far as we 
know, unprecedented. 



B. Higher Moments and Conductance Distributions 

Next, we present the results for the conductance dis- 
tributions in different regimes and provide a statistical 
description of the distribution functions, P(g). This will 
be done by focusing on (I) the largest square samples 
we have studied, i.e. 96 x 96 DNs, (2) three values of 
L/£ = 0.31, 1, and 4.09 that are typical representations 
of the metallic, crossover, and insulating regimes. 

(i) In the metallic regime, the distribution of the con- 
ductance for 5, 000 samples is shown in Fig. [?]. The solid 
line is a fit to the Gaussian distribution 



Pis) 



1 



(53) 



The central-moments, < \Sg\ n >, are computed for up 
to n — 10 from the data and compared to those deter- 
mined by the Gaussian fit in Eq. (|53|). The good agree- 
ment shows that the conductances are close to being nor- 
mal distributed in the metallic regime. The deviation at 
large values of n from the Gaussian moments observable 
in Fig. (^) is indicative of the log-normal tails in the con- 
ductance distribution function IB ill. 



< g >oc {H/Lf/ 2 e- L / 2 ^ 



(51) 




and it is therefore not meaningful to focus on the absolute 
values of the variance. In fact, a well known result is that 
the relative fluctuations are exponentially large, 



FIG. 7. The conductance distribution P(g) in the 
quasi-lD metallic regime at L/£ = 0.31. The solid line 
is a Gaussian fit to the histogram. 
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FIG. 8. Comparison of the computed central moments 
of the conductance (solid circles) with those determined 
by the Gaussian fit (solid line) in Fig. |^. 

(ii) In the insulating regime, we plot the distribution 
of the logarithm of the conductance in Fig. || obtained 
from 5,000 samples. Here the solid line is a fit to the 
log-normal distribution, 

1 (log a- (log g)) 2 

P(logg) = e 2<5iog2 s) _ (54) 

^2ir(Slog 2 g) 

The higher central-moments obtained from our data 
and the fitted log-normal distribution are compared in 
Fig. 0. 



0.15 




5 10 15 



in(g) 



follows a log-normal distribution in the quasi-lD insu- 
lating regime. Thus, instead of the conductance, it is 
the logarithm of the conductance that is a self-averaging 
quantity. 




10° 



1 10 
n 

FIG. 10. Comparison of the calculated central mo- 
ments of Ing (solid circles) with those from the fitted 
log-normal distribution (solid line) in Fig. []. 

(iii) It is important to ask how the statistical distri- 
bution changes in the crossover regime from the metal- 
lic and insulating behaviors. In this case, we focus on 
L/£ = 1 and a = 0.5. The conductances of a large 
ensemble of 10, 000 96 x 96 samples were calculated to 
achieve better statistics. The conductance distribution, 
P(g), is plotted in Fig. [n]. Remarkably, the conduc- 
tance is broadly distributed between and 1, analogous 
to the behavior of P(g) found in the Chalker-Coddington 
network model at the critical point of the integer quan- 
tum Hall transition pjjp7| , pc|. The n-th order central- 
moments are shown in Fig,^12| up to n — 10, which can 
be very well fitted by 

< Sg n >=av n e un \ (55) 
with (a, v, u) = (0.81 ± .01, 0.29 ± .01, 0.050 ± .001). 



FIG. 9. The distribution of the logarithm of the con- 
ductance in the quasi-lD insulating regime at L/£ = 4.09. 
The solid line is a fit to the log-normal distribution. 

The excellent agreement confirms that the conductance 
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FIG. 11. The conductance distribution in the crossover 
regime at L/f; = 1.0. Notice that P(g) is broadly dis- 
tributed between and 1. 
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FIG. 12. The n-th central moments of the conduc- 
tances in Fig. [n] The solid line is a fit to Eq. 

Since Eq. ( pq) is indicative of a distribution with log- 
normal tails plj , in Fig. [l3], we plot the distribution of 
the logarithm of the conductance, P[\og(g)]. The latter 
turns out to show a skewed log-normal distribution as a 
result of the sharp fall off of P(g) close to g = 1. The 
conductance can be described surprisingly well by a log- 
normal distribution as shown by the solid line in Fig. |l3j. 



FIG. 13. The distribution of ln(<;) in the crossover 
regime (see Fig. [ll]) showing the skewed log-normal dis- 
tribution (solid line). 

We point out the interesting similarity between the 
above features in the crossover regime (iii) and those 
found at the 2D integer quantum Hall transitions pl| , ^7| . 
A possible scenario to explain this is the following. The 
critical properties of the quantum Hall transition are de- 
scribed by those of the Chalker-Coddington (CC) net- 
work model at criticality |I^j20| , f32"[ . As pointed out 
in Refs [j32|j2l|| , there exists a critical manifold, corre- 
sponding to a line of fixed points, for the CC network 
model that has the same tunneling parameter at each 
node (see Eq. (||)). One can show that, as the tunnel- 
ing parameter is varied along this line, the critical con- 
ductance exhibits a crossover from 2D metallic to quasi- 
1D insulating behaviors p3| ]. The conventionally stud- 
ied quantum Hall transition corresponds to a particular 
point on the critical line selected by requiring an addi- 
tional 90 degrees rotational invariance of the network. 
The latter happens to be in the crossover regime with 
L/£ = ln(l + V5) = 0.8814 |o| and a c = 0.58 ±.03. It is 
possible that the distribution of the conductance at the 
critical point in the regime of crossover between metallic 
and insulating behaviors possesses the universal features 
that have been observed in both of the systems discussed 
above. 



C. Level Statistics 

It is well known that the difference in the statistical 
properties of the conductance in disordered metals and 
insulators originates from that of the statistical distri- 
bution of the eigenstates. A natural quantity to study 
is the distribution of the normalized (by the averaged 
value) energy level spacings, often called the level statis- 
tics. Traditionally, one applies the random matrix theory 
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developed by Wigner to study the level spacing distribu- 
tion of an ensemble of random Hamiltonian matrices [j34| . 
More recently, it has become clear, by the Coulomb gas 
analogy and the maximum entropy hypothesis P2u35[ , 
that the same theory describes the level statistics of the 
eigenvalues of the transfer matrix. The latter, in con- 
trast to those of the Hamiltonian matrix, correspond to 
the scattering eigenstates and are directly related to the 
Landauer formulation of the conductance by Eq. (21). 

Here we present our analysis of the transfer matrix 
level statistic in the quasi-lD metallic, insulating, and the 
crossover regimes of the multilayer quantum Hall surface 
states. Following Eq. (p0|), let us introduce the normal- 
ized eigenvalues for a given sample, 



~, i=l,...,C. 



(56) 



The normalized level spacing between adjacent eigenval- 
ues is given by, 



£t+i — &t 

where A is the averaged level spacing, 
C-i 



A 



C 



(57) 



(58) 



such that the distribution P(j') is uniform. The level 
spacings are then obtained by replacing 7 with 7' in 
Eq. ©. 

In Fig. yr[ the typical P(s) in the quasi-lD metallic 
region is shown for L/£ = 0.31. The solid line is a fit 
to the Wigner-surmise in Eq. @ with a = 2.02 ± .02, 
(3 = 1.88 ± 0.02, in very good agreement with the values 
expected for the unitary universality class. The small 
discrepancy of (3 from 2 may be attributed to insufficient 
statistics at small s and finite size effects. 
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According to the random matrix theory, the behavior 
of the level distribution, P(s), is determined by the level 
correlations. In the metallic phase with extended eigen- 
states, the level repulsion is strong and long-ranged, P(s) 
follows the Wigner hypothesis, 



P 



metal 



(s) = As e 



Bs° 



(59) 



with a — 2. In the unitary universality class, (3 = 2, and 
A = 32 /it 2 and B = 4/ir from normalization. On the 
other hand, in the insulating regime with exponentially 
localized states, the levels are uncorrelated, the asymp- 
totic behavior of P(s) for large s follows the Poisson dis- 
tribution, 



-PinsuO) OC e 



— ks 



s » 1, 



(60) 



where k is a positive constant. However, the behavior of 
P(s) at the critical point of metal-insulator transitions, 
or in the crossover regime between metals and insula- 
tors is an unresolved issue under active current debate 

iH§- 

From the set of the calculated eigenvalues {ji}, we 
now study the properties of the distribution P(s). Notice 
that in general, the eigenvalue distribution P( 7 ) is not 
uniform. In order to obtain a meaningful P{s), one has 
to, as is normally done |4l]j39| , p0| , unfold the spectrum of 
P(7) according to 



P( 7 )d7, 



(61) 



FIG. 14. Distribution of the level spacings in the metal- 
lic region at L/£ = 0.31. The data (+) correspond to 
L — C — 96. The solid line is a fit to the Wigner-surmise 
form in Eq. 

The level statistics becomes more intricate in the 
crossover regime. In Fig. [l5[ we show P(s) at L/£ 
1 . The most interesting feature is the coexistence of 
metallic-like and insulating- like statistics. The Wigner- 
surmise in Eq. (||), with a = 2.03 ± 0.02 and (3 = 
1.88±0.02, remains a rather good description (solid-line) 
of the data so long as s is not too large. On the other 
hand, the tail at large s is clearly described by Poisson 
statistics as shown by the dotted line. The two differ- 
ent statistics merge together around s ~ 2. This kind 
of hybrid of metallic and insulating behaviors was first 
pointed out by Shklovskii, et. al. to describe the level 
statistics at the 3D Anderson metal-insulator transition 
[^,0. These results suggest that, in the crossover be- 
tween the metallic and insulating regimes, the level cor- 
relation becomes finite-ranged, which manifests itself in 
the crossover of P(s) from the correlated Wigner-surmise 
statistics to the uncorrelated Poisson statistics at large s. 
Interestingly, such behaviors are also observed at the in- 
teger quantum Hall transitions ^7j, suggesting that the 
results of Shklovskii, et. al. may be valid for all de- 
localization transitions. Analogous to the conductance 
distribution P(g) discussed above, this is likely to be a 
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general property of the critical eigenstates. 
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FIG. 15. Distribution of the level spacings in the 
crossover region at L/£ = 1. The data (+) corresponds to 
L = C = 96. The coexistence of metallic-like (solid line) 
and Poisson statistics (dashed line) at large s is evident. 

The distribution of the level spacings in the insulating 
regime is shown in Fig. |l6| for L/£ = 6.07. As expected, 
the region of the uncorrelated Poisson tail at large s ex- 
pands as a result of the reduction of the correlation range 
of the levels. This can be seen by comparing Fig. |l5| with 
Fig. |l^ at large s. One expects P(s) to eventually follow 
the Poisson distribution deep in the insulating regime, 
i.e. for L/£ 3> 1. It is interesting to point out that 
the coexistence of metallic and insulating-like statistics 
is very robust and extends over a wide range of L /£ val- 
ues. 
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FIG. 16. Distribution of the level spacings in the in- 
sulating regime at L/£ = 6.07. The data (+) corresponds 
to L = C = 96. 

V. SUMMARY AND DISCUSSIONS 

In conclusion, we have studied the transport properties 
of the chiral edge states on the surface of a multilayer 
integer quantum Hall state in the field direction. We 
emphasized the criticality of the surface state, and that 
different behaviors of coherent transport can be reached 
depending on the path along which the thermodynamic 
limit is taken. The sample aspect ratio is an important 
quantity that enters the scaling function of the conduc- 
tance and conductance fluctuations. We presented in 
detail a new and stable algorithm for large scale con- 
ductance calculations in the transfer matrix approach to 
the Landauer formulation. This algorithm allowed, for 
the first time, a FFS analysis of the conductance and 
conductance fluctuations in the thermodynamic limit of 
the DN model of the chiral surface state. The transport 
properties, in the thermodynamic limit approached with 
a fixed aspect ratio, resemble that of a quasi-lD conduc- 
tor, showing a smooth crossover between metallic and 
insulating behaviors as a function of the interlayer tun- 
neling. The asymptotic scaling function of conductance 
and its variance are found to be in remarkable agreement 
with the analytical functions obtained using the super- 
symmetric nonlinear cr-model |l7j . The statistics of the 
two-terminal conductance is found to follow the normal 
distribution in the metallic and log-normal distribution 
in the insulating regime. In the crossover regime, the con- 
ductance is broadly distributed between zero and e 2 /h, 
which is well described by a highly skewed log-normal dis- 
tribution, similar to that found at the 2D quantum Hall 
transitions (2lJ. We also presented, for the first time, 
a detailed study of the level statistics in the eigenvalue 
spectrum of the transfer matrix. While the latter in the 
metallic regime follows the Wigner-surmise in the unitary 
universality class, coexistence of correlated metallic and 
uncorrelated Poisson statistics characteristic of an insu- 
lating state is found to describe the data in the crossover 
regime. We interpret the latter as a manifestation of 
the finite range nature of the level correlations, which 
emerges in the crossover to the insulating regime. 

Finally, we briefly discuss the transport behaviors in 
the 2D metallic regime (see Fig. ||). As discussed below 
Eq. (||), in the thermodynamic limit, the 2D chiral metal 
survives only in systems with a vanishing aspect ratio. 
The behavior of the conductance and conductance fluc- 
tuations can be obtained through FSS of a sequence of 
samples with aspect ratio A ~ 1/L a and a > 1. For 
example, one can approach the asymptotic limit with 
A = a/L, where a is a constant. Eqs ( |43| ) and ( fli]) would 
then lead to L/Lq = y/a/cr and L/£ = a/2aL. The 
scaling equations for the conductance and conductance 
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fluctuations in Eqs 
totic limit, 



and (H) become, in the asymp- 



< 9 > 



e a 



<V>=( T 



e 



2/i 



(62) 



It would be interesting to compare results obtained this 
way with the perturbative spin-wave expansion results 
derived by mapping the problem to a ID Su(n,n) quan- 
tum ferromagnetic spin chain in the limit n — > p[ . 
However, in the DN model, because the ballistic contri- 
bution dominates at large a (t 2 — > 1), in order to probe 
the asymptotic 2D regime of vanishing aspect ratio, one 
has to go to much larger system sizes |28| . The stable al- 
gorithm present here makes it feasible to carry out such 
numerical calculations. These results will be discussed in 
a future publication p3| . 
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